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Abstract. Scalar rational functions with a non-negative real part on the right half 
plane, called positive, are classical in the study of electrical networks, dissipative sys- 
tems, Nevanlinna-Pick interpolation and other areas. We here study generalized positive 
functions, i.e with a non-negative real part on the imaginary axis. These functions form 
a Convex Invertible Cone, cic in short, and we explore two partitionings of this set: (i) 
into (infinitely many non-invertible) convex cones of functions with prescribed poles and 
zeroes in the right half plane and (ii) each generalized positive function can be written 
as a sum of even and odd parts. The sets of even generalized positive and odd functions 
form subcics. 

It is well known that the Nevanlinna-Pick interpolation problem is not always solvable 
by positive functions. Unfortunately, there is no computationally simple procedure to 
carry out this interpolation in the framework of generalized positive functions. Through 
examples it is illustrated how the two above partitionings of generalized positive functions 
can be exploited to introduce simple ways to carry out the Nevanlinna-Pick interpolation. 

Finally we show that only some of these properties are carried over to rational gener- 
alized bounded functions, mapping the imaginary axis to the unit disk. 



1. Introduction 

1.1. Historical perspective. Functions which are analytic in the open right half-plane 
C+ and with a non-negative real part there 



(1.1) Rep(s)>0, s£C + , 

here denoted by V, play an important in the theory of electrical networks they were first 
studied around 1930 by W. Cauer, [22], [23] and O. Brune (who first coined the name 
positive for such functions; see [20, Definition 1, p. 25], see also [H]). These functions 
also serve as the corner stone of the theory of linear dissipative systems (a.k.a absolutely 
stable), see e.g. PH Theorem 2.7.1], [13 3.18], [S] and [17]. 
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One can weaken condition (II. ip and assume that a function p is analytic almost every- 
where on the imaginary axis satisfying 

Rep(s)>0, s G iR. 

These functions will be called generalized positive and thus denoted by QV. They were 
first addressed more than forty years ago by B.D.O. Anderson and J.B. Moore in 

We shall denote by C_ the open left half of the complex plane. We also denote the 
closed right half plane by C+ (= iRU C+). We here consider the field of scalar rational 
functions of a complex variable s with complex coefficients. Throughout this work we 
denote by QV the set of scalar rational generalized positive functions and by V its 
subset of positive functions. 

In the sequel, we shall find it convenient to denote for an arbitrary rational function g(s), 

g*(s) := GK-Or- 

The following result appeared in [8] (to help the reader we shall use p G V and tp G QV): 

Theorem 1.1. A rational function ip{s) is generalized positive if and only if it admits 
the factorization 

(1.2) ^{s)=g{s)p{s)g*(s), s£C, 

where p G V and g is rational such that both g and g^ 1 analytic in C_ and Q 
deg(g) G [0, deg(V>)]. 

Moreover, one can alwayto find s Q G iR so that in (jl.2p the functions p(s) and g(s) 
are uniquely determined by ^ ip(s ) = p(s Q ) and g(s Q ) = 1. 

Factorization results of the form of (jl.2p are well known in other frameworks. To name 
but three: 

(i) Schur functions, analytically mapping the open unit disk to its closure, rational gener- 
alized Schur functions mapping the unit circle to the closed unit disk, see e.g. [3§], [23], 
PQ. Factorization result of generalized Schur functions appeared in [381 Theorem 3.2]. 

(ii) Caratheodory functions, mapping the open unit disk to C+ and the rational gener- 
alized (=pseudo) Caratheodory functions, mapping the unit circle to C+. Factorization 
result of generalized Caratheodory functions appeared in [28, Theorem 3.1]. 

(iii) Nevanlinna functions, analytically mapping the closed upper half plane to its closure 
and the rational generalized Nevanlinna functions mapping the real axis to the closed up- 
per half plane. Factorization of generalized Nevanlinna functions appeared in parallel in 
[31] and [33] and further explored and extended to operator valued functions in [42], [43] 
and [4"4"j . 

An extended version of QV functions was introduced by M.G. Krein and H. Langer in a 
long and celebrated series of papers of which we cite only [38] , [39] (note that they studied 
functions meromorphic in the open upper half plane, or in the open disk). 

^Recall, the degree of a rational function is taken to be the maximum between the degrees of the 
numerator and the denominator. 
2 assuming t[>(s) ^ 
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Important part of the existing research on QV functions is neither confined to scalar 
functions nor to rational functions. Restricting the discussion to scalar rational functions 
mapping the imaginary axis into the right half-plane, enabled us in [8] to obtain an elemen- 
tary proof for the factorization (jl.2p . On the expense of generality, to keep the exposition 
simple, we here adhere to the case of scalar rational functions. 

1.2. The current work. Traditionally, QV functions were studied almost uniquely by 
mathematicians. They were addressed in the framework of not necessarily rational func- 
tions. In contrast, Electrical engineers have long been interested in rational positive func- 
tions (impedance of R-L-C electrical circuits). 

In this work we take a challenging task try to simultaneously address both audiences. 
Thus, depending on their background, some readers may find part of the statements 
nearly obvious and others not clear at all. Moreover, we try to address skeptical questions 
like: 

(An engineer): "Why should I be interested in the extension of positive functions to the 
QV framework?" 

(A mathematician): "Why should I be interested in scalar rational QV functions if the 
operator-valued non-rational case has already been addressed ?" 

As already pointed out above, positive functions have been a corner-stone in system theory. 
We thus believe that a sufficiently good motivation for a researcher in this field to explore 
QV functions, is to gain a proper perspective on the subset of V functions. In particular to 
understand which of the properties of V exist in the larger set QV and which are peculiar 
to positive functions. 

As a prime example we point out that rational QV functions (bounded at infinity) may 
be characterized through the Generalized Positive Real Lemma (in the positive case a.k.a. 
the Kalman-Yakubovich- Popov Lemma) see e.g. [TT], [32] for early accounts and most 
recently [9]. 

One can examine this gap between QV functions and its Hurwitz stable subset of V 
functions, from the Matrix Lyapunov Equation point of view. It has been of interest to 
identify which of the properties of the Lyapunov Equation for Hurwitz stable matrices 
are carried over to the general inertia case, see e.g. [361 Chapter 2 and Section 4.4] or 
the relation to the dimension of the controllable subspace in [40J . Resorting to the matrix 
Lyapunov equation was not just a metaphor, it is part of the above mentioned Generalized 
Positive Real Lemma, see e.g. [9], [11], [52] . 

Recall that convex sets are an essential ingredient in optimization. As an illustration, a 
typical control engineering problem would be: "Find among all stabilizing controllers the 
one which minimizes a certain property". Such problems are easier to address if the set 
of all controllers or closed loop systems is convex. Convex sets also serve as a model for 
uncertainty, e.g. the celebrated Kharitonov Theorem for checking the Hurwitz stability 
of a polytope of polynomials, see [16]. Hence, one is motivated in studying convex sets 
of rational functions which are Hurwitz stable (and then look for stable minimum phase). 
In Section [2] we identify maximal convex cones of rational functions with prescribed poles 
and zeroes in the open right (or left) half plane. 
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Using the above result, we introduce in Section [3] a partitioning of all QV functions into 
"small", yet maximal, convex cones denoted by QV g , with prescribed poles and zeroes in 
the right half plane. Each of these sets is a replica of V functions. 

Recall that a convex cone which in addition is closed under inversion is called a Convex 
Invertible Cone, cic in shorf] , see e.g. [25], [26] and [27]. It is easy to see that the set 
of QV functions is closed under positive scaling, summation and inversion, i.e. a cic and 
V is a subcic of it. More precisely, V is a maximal cic of functions which are analytic in 
C+, see [271 Proposition 4.1.1] and Proposition 12.21 below. 

In Section |4] partition each QV function to a sum of even and odd generalized positive 
functions. It turns out that the sets of even and odd generalized positive functions are 
subcics of QV. First, the set of all Odd functions (i.e. generalized lossless) which is of 
particular interest, is then studied. In Section [5] even QV functions are explored. As a 
by-product of this partitioning it is shown that it is only within the larger QV set that a 
positive function can always be written as a sum of even and odd part. 

From an applications point of view, Nevanlinna-Pick type interpolation was originally 
motivated by the design of the driving point impedance of R-L-C electrical circuits and 
subsequently by control, both restricted to positive functions. There are good reasons 
to study Nevanlinna-Pick interpolation over QV functions: 

Recall that in the positive case solution exists, if and only if, the Pick matrix associated 
with the data points is positive semidefinite. In the QV framework, this restriction is 
removed, namely for almost any set of data points a solution exists^- Note that in some 
applications, Hurwitz stability is not required e.g. when the data is not associated with 
an input-output dynamical system or when the system at hand is not passive. 

Nevanlinna-Pick interpolation problem of generalized Schur and Nevanlinna functions has 
been well addressed in the literature: For generalized Nevanlinna functions see e.g. [3], 
[ID], [5], [H], [291 Section 3] and [4J. For generalized Schur functions see e.g. 0, p3], [E], 
[19] . |30j and |41] Nonetheless from computational point of view the known procedure is 
involved. 

In each of the Sections: El [4] and (H we illustrate through examples how can one exploit the 
new structural results to simplify the Nevanlinna-Pick interpolation problem. A careful 
examination of this examples suggests directions for future research. Some of them are 
stated in Section [7J 

Recall that a function is called bounded, denoted by fa E B, (commonly the real case is 
addressed) if it analytically maps C + to the closed unit disk, see e.g. [121 Chapter 7], 
[17\ Section 6.5] and f g b £ QB is generalized bounded if it maps the imaginary axis to 
the closed unit disk, see e.g. [32] . It is known that through the Cayley transform positive 
functions may be identified with bounded functions. In Section [6J properties of generalized 
bounded functions, which do not trivially follow from this Cayley transform are explored. 

•^Strictly speaking, this means that whenever the inverse exists, it also belongs to the set, e.g. the set 
of positive semidefinite matrices is a cic. In contrast, the open upper half of C is not. 

^For generalized Schur function this follows from [30] and by appropriate Cayley transforms (of the 
functions and of the variable) this is true for QV functions as well. 
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In particular it is shown that one cannot easily mimic Section [3] to obtain a partitioning 
of all rational generalized bounded functions to a union of sets with prescribed poles and 
zeroes outside the unit disk. 

2. Maximal convex sets of rational functions with prescribed poles and 

zeroes in c + or in c_ 

In this section we consider poles and zeroes of sums of rational functions. Up to possible 
cancellations, poles of a sum are the union of the poles original functions. However, in 
general little can be said about the zeroes of a sum. We now characterize maximal convex 
sets of rational functions with prescribed poles and zeroes in C+ (or in C_). 
We begin with some preliminaries. We find it convenient to define the following sets 

^2 -|\ Q— '■= all rational functions with poles and zeroes in C_ , 

Q + := all rational functions with poles and zeroes in C+ . 

Note that g £ Q_ is equivalent to G Q + . 

Example 2.1. All degree one real functions with poles and zeroes in the C + are given by 

(2.2) Qi = \ aS ~\ : ab > 0, cd > 0, ad ^ be. 

[cs — d 

Now = | jifijpg : ab > 0, cd > 0, ad ^ be. | , and all degree one real functions within 

V is a subset of where in addition a, c > 0. Indeed, up to inversion, all degree one 
real functions in V are of the form as + b or - + b with a, b > (in electrical circuits 
terminology, the driving point impedance of R-L or R-C networks, respectively). □ 

Next, for a given g + in Q + , see (|2.ip . let Q g+ be the set of all rational functions with 
prescribed poles and zeroes in C+ given by, 

~ I TL [Si 

(2.3) Q g+ ■— \ , (\ 9+( s ) '■ n-{s),d-(s) polynomials with roots in C_ 

I d— ( s ) 

Obviously, if in (|2.3[) g G Q g+ satisfies = const, then in fact g G G+ (12. ip . The 

set Q g+ can not be convex as both — g+(s) and sg + (s) belong to it, but their sum has 
an additional zero in C+. Yet, it is of interest to identify maximal convex subsets of Q g+ 
in (|2.3p . To this end, denote by Q Q the set of rational functions, which along with their 
inverses, are analytic in both open half planes. Namely, whenever g Q G Q , it is of the 
form 



(.s - /rj)' IJ 



(2-4) g ( S ) = cH(t 

i=i 

where rj G M are all distinct, rjj are integers (positive or negative) and c G C. We shall 

o 

use the convention that J} = 1, so that also g(s) = const, belongs to this set. 

i 

For three given functions g + G G+ , g~ G Q- see (|2.1[) . and g Q see (12. 41) . define, 
(2-5) Q g+ , g _, go := {g+(s)p(s)g (s)g-(s) : p G V }. 
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By construction, in C + the poles and zeroes of all functional in G g+ , g _ t go are exactly 
those of g+(s). On iR poles and zeroes are almost fixed in the following sense. Consid- 
ering dZ3]), functions in G g+ , g_, g have factors of the form (s — ir,-)^ + , (s — irj) 11 ^ 1 
or (s — irj) r>i depending on p(s) having at s = irj, a zero, a pole, neither zero nor pole, 
respectively. 

We can now describe convex sets of functions with prescribed poles and zeroes in C+ 
and on iK, almost prescribed in the above sense. 

Proposition 2.2. The following is true: 

(i) The set V is a maximal convex invertible cone, cic, of rational functions analytic 
in C+. 

(ii) G g+ ,g_,g in (|2.5p . is a maximal convex set of rational functions with prescribed 
poles and zeroes in C+. 

In fact, if cj) Gg + ,g-,g a , one can always find tp G G g+ , g _, go so that (<p + ip) g" Q g+ 

•CEO). 

The fact that the set V is a convex invertible cone, cic, is classical, see e.g. |17l 5.6]. The 
result in item (i) is a small variation of \27\ Proposition 4.1.1]. 

Proof : (i) We first show that if h(s) is a non-positive function, one can always find 
a positive function p so that h + p has a zero in C+. Indeed, let h $ V be given. 
By definition there are points in C+ which are mapped by h(s) to C_, i.e. there exist 
a,7>0, (3,5 £R so that h(s)\ a=a+i0 = -7 + iS. Take now p(s) = %s - + 5). Then 
clearly p G V and (p + /i) (s) [ a=a+j|S = 0, i.e. a zero in C + . Next, note that is not 

analytic in C+. Since V is closed under inversion (i.e. p 6 P is equivalent to = G "P), 
this part is established. 

(ii) Let g- £ G- , g+ G ^+ see (12. ip and g„ G £7 , see (I2.4p . be given and let <f>{s) be a 
rational function not in G g+) g 0) g_. To avoid triviality assume that <f> G {Gg + \ G g+ , g _, go }- 
Next, denote, /i(s) = 5o(s)~ 1 5+(s)"V(s)5 r -(s) _1 - Then by (J23|), /i P (else ^> would 
have been in G g+ ,g_, go )- Take now tf)(s) = g + (s)g (s)p(s)g-(s) with the above g+(s), 
g (s), g~(s) and p(s) as in part (a) of this proof. By construction, p G V and thus 
^ £ G g+ ,g_,g , but ((f) + ip) since this function has an additional zero in C+. (If this 
additional zero coincides with an existing pole, both in C+, still (cj) + -ip) g" ^ s+ ). Thus, 
this part of the claim is established and the proof is complete. □ 

As an illustration we have the following. 

Example 2.3. Take g+(s), g (s) and g~(s) in (|2.5p to be fixed polynomials. A maximal 
(up to scaling) convex set of polynomials whose roots in C+ are those of 5+, is given by 

{g+(s)(s + a)g (s)g-(s) : a G C + }. 

Indeed, p(s) = s + a, a G C + are the only polynomials in V. □ 



'excluding the zero function 
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We conclude this section by stating the analogous results for the left half plane. First, we 
denote by "P* the set of para-positive functions, 

(2.6) V* := { ip : ip* eV}. 

Thus, functions in V^ map C_ to C+. In particular, V^ C QV. We can now state 
results which are dual to Proposition 12. 21 

Observation 2.4. The following is true: 

(i) The set V& is a maximal convex invertible cone, cic, of rational functions analytic 
in C_. 

(ii) Let g_ G g + G Q+ and g a G Q Q be given. The set 

{g + (s)p#(s)g (s)g_(s) : p € V } 

is a maximal convex set whose poles and zeroes in C_ are precisely those of 
9-(s). 

3. Convex partitioning of QV functions 

We now address ourselves to subsets of generalized positive functions within Q g+: g _ t go 
in (|2.5|) . namely sets of the form Q g+j g> 9o n QV. To this end, we introduce the following 
set, using (pT|) and ((23]), 

(3.1) Q+ ■= {g + {s)g (s) : g + G Q+, g Q G Q Q }. 

Note that g G Q+ means that both g and g^ 1 are analytic in C_. For example, all 
degree one functions in Q + are given by (the real subset Q\ was described in (|2.2p ). 

(3.2) Q = l^zl : Re( a *b) > 0, Re(c*d) > 0, ad ^ be] . 

[cs — d J 

Using this notation, we shall hereafter simply write 

<3V 9 ■ 0:,.., OV. 
By Theorem ll.il for given g G Q+ this can be equivalently written as 

(3.3) QV g = { gpg* : p£V}. 

For a given g £ Q + , the set C/'Pg is a replica of P. Nevertheless, the picture in QV g is 
richer. 



Example 3.1. We here illustrate two properties of the set QV g (j3.3|) . where g G Q+ is 
given: 

(a) If ipj = gpjg^ with g fixed, degpi > degp2 ; does not always imply deg -01 > deg if)2 ■ 

(b) In this set, it is only in C+ that the poles and zeroes are fixed. On iM they are 
almost prescribed, and in C_ they are not fixed. 

Take g(s) = ^ {g G Qi, see ([22D). Thus, QV g = { ^p(s) : pG7?}. We here 
mention, but five interesting samples, 
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p(s) ip(s) = g(s)p{s)g#{s) 

s s-2 
s+2 s 
s s-2 

(s+2) 2 s(s+2) 

1 ^ 

s+2 (s+2) 2 {s-2) 

s(s+2i) (s 2 -4)(s+2i) 
s+i s(s+i) 

(a) These five functions are ordered so that the degree of ipj(s) is non-decreasing. In 
contrast, the degree of the corresponding Pj(s) "fluctuates". 

(b) In C+, there is always a zero with a unit multiplicity at s = +2. 

On iM., there is a pole at the origin. Its generic multiplicity is two, but it may also be one 
or three (i.e. at the origin p has a pole e.g. (iv), no pole nor zero e.g. (hi), or a zero e.g. 
(i), (ii), (v) respectively), ifi may have additional poles or zeroes, see e.g. (v). 
In C_ poles and zeroes are not fixed, see e.g. the point s = —2. □ 

Using the notation of (|3.ip Theorem 1 1 . 1 1 may be formulated as saying that having ifi G QV 
is equivalent to t/)(s) = g(s)p(s)g#(s) for some g G G+ and some p G V. Thus, we can 
use the last result to introduce a convex partitioning of all QV functions. The proof is 
immediate and thus omitted. 

Observation 3.2. Let QV g be as in (j3.3[) . Then, 

(i) QV g is a convex sub-cone of QV . 

(ii) For 9l ,g 2 eg^ QV rj2 = (f)QV gi (ff. 

(iii) Let g±,g 2 G Q+ be so that g 2 ^ cgi, for some constant c, then, QV gi n QV g2 = {0}. 

(iv) QV = \J_QV g . 

(v) {QV g )- x = QV {a * ) -x. 

Proof. Items (i), (ii) and (v) are immediate from (|3.3p . Item (iv) follows from Theorem 

O along with (03]). 

As to item (iii), assume that there exits V within QV gi n QV g2 for some gi,g 2 G Q + . 
We shall find it convenient to factorize gj = g Q jg + j with j = 1,2. where g ,i,g ,2 £ Qoi 
see (|2.4p . and g+,i,g+ t2 G Q+, see (|2.ip . As poles and zeroes of ip in C + are uniquely 
determined by g + , without loss of generality one can write g\ = g ,ig+ and g 2 = g ,2g+ 
for some g + G Q+. 

Next assume that for j = 1,2 and some r G K, go,j(s) have factors (s — ir) m :> and the 
corresponding Pj(s) have factors (s — ir) i , where rrij,lj are integers (not necessarily 
positive). Then in ipj(s) the respective factors are 

(s - ir) m ^(s - irfi ((s - ir) m i)* = (-l) m ^(s - irf m ^ l K 

This implies that: (i) 2m\ + 1% = 2m 2 + Z2 and (ii) mi — m-2 = 2k for some integer k. 
Namely, l 2 — h = 2 (mi — ma) = 4fc. Now recall that on the imaginary axis poles and 
zeroes of positive functions are simple, see e.g. [U Theorem 2.2], i.e. 1 > This implies 
that 2 > \l 2 — But, mi 7^ m2 implies that \l 2 — h\ > 4. Hence, one can now conclude 



(i) 
(ii) 

(iii) 
(iv) 
(v) 
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that mi = rri2, l\ = l 2 and since r was arbitrary (up to a constant) 50,1 = 5o,2j which in 
turn implies g\ = g 2 (up to a constant) and the proof is complete. □ 

This partitioning of QV functions rightfully seems straightforward. In contrast, at the 
end of Section [6j we show that partitioning of QB, generalized bounded functions (or 
generalized Schur functions) in the spirit of Observation 13.21 can not be easily mimicked. 

For a given g G Q+ we now wish to identify minimal degree functions within QV g . 

Proposition 3.3. The following is true. 

(i) g G Q + can always be factored as g(s) = c Y[ 4>j(s) with ^{s) positive, see 

i=i 

(l2i)D . and c G C. 

(ii) Among all possible factorizations of g G Q+ let us choose g(s) = gi(s)g2(s) so 
that gY(s) is positive and deg(gi) is maximal. 

Then, ggl[ is the minimal degree function in QV g . 

Proof : Item (i) is immediate from item (i) in Observation 12.41 Specifically, g G Q+ can 

n(«-*i) 

always be written as g = frr^^j with 2^,7^ G C+. Note now that (s — Zj)# = —(s+Zj) 

k 

and r^z^w = s+n* ' Namely, up to sign change, this 5* is a product of degree one 
positive functions. 

Item (ii) stems from item (i) and (13.3P noting that ip G QV g can always be written as, 

"000 = g{s)p{s)g*{s) = gi(s)g 2 (s)p(s)gf{s)gf(s), 

so that gf(s) is positive, see e.g. Example 12. 1[ Thus, to reduce the degree of the above 
ib(s) choose p(s) = 2, ; so that 

g*(s) 

i>(s) = 9i{s)g 2 (s)p(s)gY(s)g^ (s)i _ x = gi{s)g 2 {s)gJ[ (s) = g{s)g?[ (s) 

of 

and the construction is complete. □ 
The above construction is illustrated in part (a) of Example 13.11 

Within the set V, the Nevanlinna-Pick interpolation problem is classical, see e.g. |144 
Chapter 18], [38J. Within QV, variants of this interpolation problem are well studied, 
see e.g. [3], [TO], [5], [18], [29], Section 3] and [4] for generalized Nevanlinna functions and 
for generalized Schur functions see e.g. [TO], [TO], [TO], and [3D] and [2]. Nonetheless from 
computational point of view the procedure is involved. As an intermediate step, in the 
following example we illustrate the fact that within the set QV g , namely when g G Q+ is 
fixed, the Nevannlin-Pick interpolation problem reduces to the classical version within the 
set V, which is computationally well established. 

Example 3.4. We here illustrate a Nevanlinna-Pick interpolation scheme within the set 
QV g . We look for ip G QV g so that 

i){s)\ B=1 = l and V(s)| s=2 =4. 
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(As the associated Pick matrix is ^5 |J, its determinant is negative, so there is no if) E V). 

Tak^| g(s) = 7 ^ 3s i.e. a right half plane pole at s = | and a zero at infinity. Thus, 
g(s)g*(s) = . First denote 

wi := g{s)g#(s) u ^ 1 =§ and w 2 := g(s)g# (s)\ s=2 = if . 
Exploiting the C/T'g structure, see (|3.3p . we actually seek p £ V so that 

K s )| s=1 = ^ =| and p(s)| s=2 = 4 = f • 
But this is a classical Nevanlinna-Pick problem and the resulting Pick matrix is, IT = ( % f§ 

V 12 s 

As this II is positive definite, there are infinitely many solutions (using II, they can all 
be parameterized, see e.g. [HI Chapter 18]). 

Take for instance two degree three interpolating functions, p a (s) = ^fejijj an d 
Pb(s) = I + 3( s %+2) ) • T ne resulting interpolating functions in QV g are 

^o(a) = 5(s)Pa(s)5(s) # = ( a 42)(49"-L M ) ' 
V»6(s) = 5 , (s)Pfc('S)5(s) # = ( ^ 2 + 2) {49_t2) • 

Both if) a (s) and ^(s) are of degree four. 

Recall that the set of interpolating functions is convex. Take for example 

Ms) :=# (s)+ iM s ) = 9(s) = , 
to obtain a unity degree interpolating function within QV g . From item (ii) in Proposition 
13.31 it follows that in fact this is the minimal degree function within QV g and in particular 
the minimal degree interpolating function. □ 

4. Odd functions - a suBcic of generalized positive functions 

As already mentioned, it is easy to see that the set of QV functions is closed under positive 
scaling, summation and inversion, i.e. a cic. In the previous section we introduced a 
partitioning of this cic into (infinitely many non invertible) convex cones of the form QV g . 
We now explore a partitioning of each generalized positive function into even and odd 
parts. It turns out that the sets of even and odd generalized positive functions are two 
subcics of QV. 

Abusing the real case terminology, for a given rational function /(s) we shall define the 
even and odd parts as 

(4.1) /even( S ) := \ (/(*) + /#(«)) /oddto == £(/(*) " / # 0»))- 

Then, we also define the sets of all even and all odd functions, 

(4.2) Even := { f(s) : f = f CVCD } Odd := { f(s) : / = / odd }. 

The following observations are almost obvious, they are stated for a comparison in the 
sequel. 

Proposition 4.1. Let the sets Even and Odd be as in (|4.ip and (|4.2|) . 



•'Recall that we assume g is prescribed. 
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(i) Even and Odd are convex invertible cones, cics of rational functions. 

(ii) Let f, g be rational functions. If g G Even then, 

(/sOeven = /even5 (/#)odd = foddd- 

Conversely, if (fg) eve n = fewenQ and f ^ 0, then g G Even. 

Proof (i) Indeed, positive scaling and summation are obvious. As to inversion note 
that if /(s) = -^fy, with n, d polynomials, then / G Even is equivalent to ndfi = n#d, 
which in turn means, / _1 G Even. The reasoning for / G Odd is similar and thus 
omitted. 

(ii) A straightforward calculation shows that having (fg) ev en = fevend is equivalent to 
f#g = fftgft, which in turn for / ^ means g G Even. □ 

From item (i) in Proposition 12.21 and item (i) in Proposition 14.11 it respectively follows 
that V and Odd are cics. Recall that a non-empty intersection of cics is a cic, see 
[27\ Observation 2.1]. Thus, of particular interest is the subcic of all positive odd rational 
functions VO := VnOdd. The real subset of "PC functions are sometimes called "lossless", 
or "Foster" and they correspond to L-C circuits, see e.g. [12] . [T7] , 

VO functions can be parameterized, see e.g. p2J 5.13], as 



(4.3) VO :-- 



p(s) = ir Q + a Q s + — —. — : a Q > 0, a 7 - > 0, r, G 



3>1 S 



Combining Theorem 11.11 together with (|4.3|) we have the following. 

Proposition 4.2. Let the set of odd functions, Odd, be as in (|4.1|) and (|4.2|) . The 

following statements are true. 

(i) Odd = GVr\Odd. 

(ii) tp G Odd if and only if tfj maps the imaginary axis to itself. 
(hi) The set Odd is closed (excluding the zero function) under: 

(a) real scaling, (b) addition, (c) inversion (d) composition and (e) the product of 

2m+l 

an odd number of elements, i.e. \\ ipj(s) with ipj G Odd, m = 0, 1, ... 

3 

(iv) t/j G Odd can always be written as 

a ; 



i)(s) = g(s) ir Q + a Q s + ^ J '. g*{s) 



with a a > 0, aj > 0, rj G K and g G G+ (|3.1j) . 
Proof : (i), (ii) Recall that for an arbitrary function /, 

(4-4) (f*(s)) i =(f(s) u J*. 



Thus, whenever / G Odd it maps iM. to itself and hence it is a QV function. Next, 
from (14. ip and (|4.4|) it follows that / eve n! the even part of an arbitrary /, maps iR to 
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R. Thus, if / = /even + /odd maps zR to £R, it follows that / eV en(s)i ^ = 0, which in 
turn means that feven(s) = 0. 
Item (hi) follows from item (ii). 

Item (iv) follows from item (i) together with (I4.3P and Theorem ll.il □ 
We now explore the structure of Odd from a geometric point of view. 

Observation 4.3. For a given set of rational functions G (G ^ {0}), denote by F and 
H the following sets, 

F := {g- 1 : g e G} H := { h = g 2 : g G G} 

G C Odd if and only if 

(4.5) Re^f(s)h(s)) UeiM )=0 V/eF Vkff. 

Proof : The claim relies on Proposition 14.21 item (ii). First, if g a (s) and <?&(s) are odd 
functions then Im {g 2 a (s) u J = and Re {g^(s) u J = 0. Thus, Re ([gl(s) g ; 1 (s)) u J = 0, 
i.e. (|4.5p is satisfied. 

Conversely, if (14. 5p is satisfied for all f £ F and all h £ if, it in particular holds for 
/ = <7 _1 and h = g 2 with the same 5. This implies that fh = g, which means g G Ocid □ 

Note that for a pair of functions /, h one can define a function- valued inner product 
< f,h> := Re (f(s)h # (s)) (in the sense that < rf, h >= r < /, h > for r G R). Now, 
as (/(s)/i(s))| sg . R = (/(s)/i # (s))| , equation (|4.5j) can be written as < f,h >| s6iR = 0. 
Namely, one can interpret Observation 14.31 as saying that 0^g£ Odd is equivalent to 
having the restriction to the imaginary axis of g and g 2 , orthogonal in the above inner 
product. 

Based on Observation 13.21 and Proposition 14.21 we can now introduce a convex partitioning 
of all Odd functions. For a fixed g G G+ the set QV g was defined in (|3.3p . We now 
consider the set odd functions in it: 

Oddg := QV g n Odd. 

This subset is given by 

(4.6) Oddg = [j g lir + a s + Y^ 

a >Q,a,j>Q,rjEM \ j>l 

Observation 4.4. Let Oddg be as in (|4.6p . Then, 

Odd = [j Oddg . 

Recall that in [3S] it was shown that (up to possibly compromising the minimal degree 
interpolating function) without loss of generality the classical Nevanlinna-Pick interpola- 
tion may be confined to V Pi Odd. Similarly, for a fixed g G G+ interpolation of QVg 
functions may be confined to Odd g functions. This is illustrated next. 
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Example 4.5. In Example 13.41 we studied a variant of Nevanlinna-Pick interpolation 
problem within the set QV ' g with g(s) = 7 ^ s , and looked for a function ip so that 

i>(s)\ s=1 = 1 and V>0)u =2 = 4 - 

It turned out that this was equivalent for a classical Nevanlinna-Pick interpolation problem 
of searching p € V so that 

P{s)\ s=1 =f and p(s)| s=2 =f • 
Following the above analysis, looking for ip S Odd g is equivalent to restricting the classical 
Nevanlinna-Pick search for p£?fl Odd. From [48J it follows that this restriction, does 
not limit the solvability of the problem. In fact, there are still infinitely many solutions. 
We here mention two solutions p a (s) = 3 l^2+2) (from Example l3.4p andpd(s) = |(8s + ~). 
The resulting interpolating functions in Odd g are 

^o(fl) = g(s)Pa{s)g*(s) = ( 8 *+2)(49-98*) > 
Ms) = g(s)Pd(s)g*(s) = 3a(49_92) • 

The function ip a (s) is of degree four while ipd(s) is of degree three. 

Following item (iii)(a) in Proposition l4.2l taking rGR as parameter, rip a (s) + (l—r)ip c i(s) 
forms a variety of interpolating functions within the same Odd g . 

Finally, a straightforward use of (j4.3|) reveals that the above ipd is a minimal degree 
interpolation function in Odd g . □ 

5. Even generalized positive functions 

In Proposition 14.21 we showed that all odd functions are generalized positive functions and 

characterized them. We now characterize QVE := QV Pi Even, the subset of even functions 

o 

within QV. We shall use again the convention that F| = 1. 

l 

Proposition 5.1. The following are equivalent 

(i) i/j £ QV£. 

m 

n (i-mi+cs-^-) 2 )) 

(ii) ip(s) = c- 3 -^ with c> 0, aj,j k e(Q, 1], (3j,5 k £R. 

n (i-7fc(i+(s-^) 2 )) 
*;=i 

(iii) ip(s) = g(s)g#(s) for some rational g(s). 

(iv) ip(s) maps iM. to M+. 

(v) i/j £ QV maps iR to R. 

Proof Any even function maps iM. to R, thus (i) implies (v). 

(v) <^=^ (iv). From (v) it follows that ip maps iR to C+ n R, thus in fact to R + , 
so this part is established. 

(iv) <J=^ (Hi). From Theorem 11.11 together with with fact that 

g(s)^(s)g*(s)\ si _ m = g(s)ip(s) (g(s))*, it follows that ip(s) = g(s)p(s)g#(s), where p G V 
maps iR to R + . Now, up to non-negative scaling, p(s) = 1 in whole C, is the only 
function which achieves that. 
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m 

n(«-*j) 

(m) (ii). Denote, g(s) = with c,zj,pk G C. Thus, 

n («-p*) 



5 ( S )<?(s) # = |cl 2 • ^ 



n -s 2 + s (^-z*) + i^i 



n 

n (-* 2 + *(p*-i$) + bfcl 2 ) 

fc=l 
m 

n ((Refe)) 2 - ( S - ilmfe)) 2 ) 

i2 

n 

II ((Re( K )) 2 -( S -iIm(p fc ))2) 
fc=l 

m 

n (1 + (Refe)) 2 - (1 + (s - Hm( Zj )) 2 )) 

1 2 £^ 

n 

II (1 + (Re(p fc )) 2 - (1 + (s - ilm(p k )) 2 )) 
fc=i 

n (i + (Re(^)) 2 ) n (i - T+(Rfe)F (1 + (s - iIm ^)) 2 ) 

1 2 i =1 i =1 



n (i + (Re(p fc )) 2 ) n (i - OTfeiF^ 1 + ( s - iIm (^)) 2 : 

m 

II (l+(Re(^))2) 

Denoting c := \c\ 2 ■ ^— — — , ay = i+^t^ , & = M*i), 7fc = i+(Re( Pfc ))^ and 
fc=i 

5fe = Im(pfc) completes the construction. 

As trivially (ii) ==>• (i), the claim is established. □ 

From Proposition 15.11 it follows that ip(s) is in QVE may be characterized as tp = gg& . 
Note however that this factorization is non-unique, namely, one can have g\ ^ gi and still 
gigf" = g2df ■ A characterization of all these factorizations is given in |34j . 
One can now state several properties of QVE, the subset of even functions within QV. 

Proposition 5.2. Let ip be a rational function. The following statements are true. 

(i) tpegr v even G QV. 

(ii) QVE is a subcic of QV. 

(hi) QV n Even is a multiplicative group. 

(iv) g G QVE <=^> g ■ QV C QV. 

jj // 

(v) For arbitrary QVE functions g\{s)g{{s) , ... , g m {s)gm(s), there always exists 
gg# g t/'P£ so £/ia£ 

m 

Moreover, one can always take g £ Q+. 
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(vi) Let tp be the composition function ip(s) ■ = p(g(s)) where g E QVS and p E V. 
Then, tp E QV. If in addition p leaves the real axis invariant (e.g. p is real), 
then V E QVS . 

Proof (i) This follows from the fact that tp E QV <=^ tp* E QV. 

(ii) Recall that the set QV is a Convex Invertible Cone, the claim is immediate from 
Proposition I5.1f iv). Alternatively, QVS is a non-empty intersection of two cics and thus 
a subcic, [STJ Observation 2.1]. 

(iii) Is immediate from items (ii) or (iv) in Proposition 15.11 

(iv) If / is a generalized positive function, from item (i) we know that /even 6 QV. Now if 
g E QVS from item (ii) it follows that (/even*?) E QVS. Next, from Proposition I4.1f ii) it 
follows that (fg)even = /even9 and hence, (fg)even E QVS. Using again item (i) implies, 
(fg) E QV n Sven. 

For the other direction all we need to show is that if g is a non-even function within 
QV, one can always find / E QV so that (fg) QV. Indeed assume that g is so that 
g(s)\ =iu = a + ib with a > and / 6 G 1 for some u Q E K. Then, taking the (odd) 
QV function / = b(uj s + i) reveals that Re ( y f{s)g( y s)\ _. ) = — 6 2 (1 + w 2 ) and thus 

gf^QV. 

(v) Is immediate from item (ii) here together with item (iii) in Proposition 15. 1[ 

(vi) By construction g maps iM to M+ and in turn p maps M+ to C+, thus, ip 
maps iM. to C+. If in addition p maps K.+ to R + , ip maps iM. to M+, so the claim 
is established. □ 

The convexity of the set QV and of its subset of QVS functions, see item (v) in Proposition 
15.21 may be exploited to introduce a straightforward scheme of solving the Nevanlinna-Pick 
interpolation problem. 

Example 5.3. For simplicity we consider the real case, 
a. We first look for a real polynomial gg& e QVS so that 

(5.1) g(s)g(s)# u=±i = 1, g(s)g(s)# u=±2 = 4, 

(Note that f(s) = s 2 is an interpolating polynomial, but only — / E QVS). 
Consider the following real QVS polynomials, 

gi (s)g#(s) = (A-s 2 )(9-s 2 M -s 2 ) 
g2(s)g*(s) = (l-s 2 )(9-s 2 )± (f-s 2 ) 
g 3 (s)gf(s) = (l-s 2 )(4-s 2 )^ (90 - S 2 ) 

It is easy to verify that, 

s = ±1 ±2 ±3 

gigf =i o o 

g 2 g# =0 4 

g 3 gf =0 9 

Thus, using item (v) in Proposition [521 55* = 9i9\ + "5252 + dzof ls a rea ^ 0V£ polynomial 
satisfying (15. lj) . 



5(s)5(*) V ±3 = 9- 
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b. Using part a, we now look for a real polynomial ip G QV so that 

(5.2) ^s) u=1 = l, i>(s) u=2 =4, ^( S ) u=3 =9, 

(with no constraints on ip(s)\ s= _ 1 _ 2 _ 3 ). For j = 1, 2, 3 we now construct real QV polyno- 

mials of the form ipj = gjPjgJ with gjgj from part a and pj E V are of the form a J +s 

with aj > is so that one of the roots of gjgj is canceled. Indeed take, 

Pl(s) = 4g , P2(S) = Jjf , P3 (S) = §£2±f . 

Thus, one obtains, 

Vi(*) = (4- S 2 )(9- S 2 )(^ 
Ms) = (1 - s 2 )(9 - s 2 )|(^+l)(^ -s) 
Ms) = (1 - s 2 )(4 - s 2 )^ (VW + I)(3v/10 - s) 
It is easy to verify that, 

s = 1 2 3 

V>i = i o o 

V> 2 = 4 
V> 3 = 9. 

Thus, taking ip = ipi + ip2 + ^3 satisfies (|5.2p . 

Roughly speaking, the simplicity of this scheme of constructing interpolating functions, 
comes on the expense of high degree. □ 

As already mentioned, the fact that the set V is a convex invertible cone cic, is classical, 
see e.g. \17\ 5.6] and item (i) in Proposition 12.21 Recall that real V functions are 
identified with the driving point impedance of R-L-C electrical circuits [12], [T7], |20j . 
|21j . In fact, in the framework of R-L-C electrical circuits the three cic operations of 
positive scaling, summation and inversion have the physical interpretation of transformer 
ratio, series connection of impedances and impedance/admittance duality, respectively. 
Moreover, recall that VnSven is associated with resistive circuits and VO with reactive 
networks (L-C). However, not every network can be realized as a series connection of a 
resistive and a reactive circuits. Namely, it is only over QV that the partitioning of a 
positive function into even and odd parts is always possible. 



Recall that in contrast to (14. 3D . the set V n Even is almost empty, i.e. up to positive 
scaling it consists of a single function, p(s) = 1. We now show that if one is interested 
in even-odd partitioning of functions, the set QV is closed, while its subset of positive 
functions is not. Namely, p e ven, the even part of a positive function p, is either a non- 
negative constant or not a positive function. One can only guarantee that p e ven 6 QVE. 
This is illustrated next. 

Example 5.4. Consider the positive (real) function: V , ( s ) = itjr; defined in the whole C. 
Then p even (s) = and fWd(s) = are so that p even G QV£ and p odd G QVnOdd, 
but neither p ev en nor p Q dd are positive. □ 

We conclude this section by introducing yet another factorization of QV functions through 
odd functions. 
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Observation 5.5. ip G QV if and only if there exist f,g£ Odd so that ip even = —f 2 
and Vodd = 9- 

Proof : Since ipodd G Odd, we only need to show that tp G QV if and only if ipeven = —f 2 
for some / G Odd. Now from item (i) in Proposition 15.21 it follows that having tp G is 
equivalent to tp evea G QV. From item (iv) in Proposition 15.11 this in turn is equivalent to 
VWen mapping iM. to K + . Using item (ii) from Proposition 14.21 completes the proof. □ 



6. Generalized bounded functions 

Recall that a function fb(s) is called bounded, denoted by /& G B, (commonly the real 
case is addressed) if it analytically maps C + to the closed unit disk, see e.g. [121 Chapter 
7], [13 Section 6.5] and f g b(s) is generalized bounded f g b G QB if it maps iR to the 
closed unit disk, see e.g. [32J. It is well known that through the Cayley transform one can 
identify positive functions with bounded functions, namely 

Nevertheless, we here focus on the less obvious analogies. In Proposition 16. II and Corollary 
below we introduce two representations of QB functions. 



Proposition 6.1. A rational f g b(s) is a generalized bounded function if and only if it is 
of the form f g b(s) = fb(s)/(3(s), where fb^B and where (3{s) is a finite Blaschke product. 

Proof: We first note that since f g b is bounded on the imaginary axis, all its singularities 
there are removable. Let w\, . . . , wg be the poles of f g b in C+, and consider the function 



Ms) = f gb ( S ) n 



s 



s + w* 

3 =\ 3 

The function f g b is analytic and bounded by 1 in modulus in C+, as is seen for example 
by the maximum modulus principle, or by direct inspection. We thus have the result with 



m = n 



3=1 J 

The converse is clear. □ 

One can characterize generalized bounded functions through the associated kernel. 
Corollary 6.2. A rational f(s) is a generalized bounded function if and only if the kernel 

WOO/M* 



kf(s,w) 



s + w* 



has a finite number of negative squares in C + \ {wi, . . . , ,W(}, where the Wj denote the 
poles of f in C + . 
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Proof: Assume that / = f g b = fb/fi- As proved in a more general context in [7\ Theorem 
6.6, p. 132], one direction follows from the equality 

1 1 
see for instance the formula on top of page 134 in [7J. 

The converse is just a particular case of the above mentioned result of Krein Langer [38, 
Theorem 3.2]. A direct proof for the rational case can also be given, but will be omitted 
here. □ 



We now turn to another representation of f g b £ QB. Obviously, (generalized) bounded 
functions and (generalized) positive functions are related through the Cayley transform 
(|6.ip . We now introduce an adapted version of this characterization. To this end, recall 
(Proposition 15. 2| ) that the set QV H Even is characterized by functions of the form gg# . 
From the above discussion one has the following. 

Corollary 6.3. A rational function f g b{s) is generalized bounded if and only if it admits 
a representation, 

(6.2) f gb (s) = (g(s)g*(s)- P (s)) (g(s)g*(s) + p( S )) _1 , 
for some p 6 V and some g(s) £ G+, (|3.ip 

Proof Indeed, from (jl.2p and (|6.ip it follows that f g b £ QB can be written as, 
f gb = {l-^){l + ^y l = {l-gpg*){l+gpg#Y l 

= (g {{g*g)- 1 - p) g*) (g {{g*g)- 1 + p) g*) ~' 
{(g*gr 1 -p){(g*g)- 1 +p) 

Now, from Proposition 15.11 and Proposition [52] it follows that ij) G GVDSven is equivalent 
to ip(s) = [g^(s)g(s)) , so up to inversion, the claim is established. □ 

It is interesting to compare Corollary 16.31 with f g b in (|6.ip 

We conclude by pointing out that there is a structural difference between QV and QB 
functions. One may be tempted to try to mimic, in the framework of QB functions, 
the convex partitioning of sets of functions of the form QB g in the spirit of (j3.3[) and 
Observations 13.21 and 14.41 se ^s where in C+ the poles and zeroes are fixed. However, 
unfortunately this is no longer true. This also prevents us from mimicking the interpolation 
over QV g to QB g functions. 

Indeed, for a given g in (|6.2p define the set QB function, 

(6.3) QB g :=\^g{s)g#{s)-p{s))[g{s)g#{s)+p{s)y 1 : P ev\. 

In the following example we show that in contrast to the set QV g in (|3.3p . within the set 
QB g neither the poles nor the zeroes in C+ are fixed. 
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I 

Example 6.4. Fix in (|6.3j) g € Q+ , see (|2.ip . namely, g(s) = § with n = n (s - 2j) 

and d= n (s — irk), where z%, ... , z\ and 7rx, ... , ir q are prescribed (not necessarily 
fe=l 

i „ <? 

distinct) points in C+. Let now, n = \\ (s — Zj — €j) and d = \\ (s — iTk — Sf.) with 

j=l k=l 

l q 

ej > 0, 6 k > 0, 1 » Z *j > °> 1 » E ^ > 0. Let pa := | and p 2 := f By 

j=i fc=l 

construction, each of the functions pi,P2 analytically maps C+ to a neighborhood of 
the point +1 and hence both are positive. Thus both 

fgb,l ■= [g(.s)g*(s)-pi(s))(g(s)g*(s)+p 1 (s)) 1 
fgb,2 ■= (g{s)g*(s) - p 2 (s)) (g(s)g#{s) +p 2 {s)y 

are in QB g . However, in C+ they do not share the poles nor the zeroes. Indeed, 
substitution yields f g ^i = (ng# — d)(ng# + d)~ l and f g \ ) ^ = {ng^ z — d){ng^ z + d)~ l . □ 

It should be pointed out that the above discussion reflects a property of QB g functions, 
independent of the choice of the representation. Indeed, if the set in (16. 3p is substituted 
by the analogous one, based on Proposition 16. 1( a conclusion, similar to that of the above 
example, is reached. 

It should be emphasized that rational generalized Nevanlinna functions, mapping the real 
axis to the upper half plane, admit a partitioning along the lines of Section [3l In contrast, 
rational generalized Schur functions mapping the unit circle to the unit disk, share the 
same difficulty as QB functions. This suggests that the complicated known scheme for 
solving Nevanlinna- Pick interpolation problem for generalized Schur functions, see e.g. 
[13], [15] and [30] and for the single point with derivatives version, [2], can not be simplified 
along the lines suggested in Example 13.41 

7. Future research 

In this work, part of ongoing research on QV functions, we concentrated on exploring 
structural properties this set. This opens the door for studying various questions and we 
here mention sample of those. First, in the framework of Nevanlinna-Pick interpolation 
problem of scalar rational QV functions. 

• Explore the question of minimal degree interpolating functions. 

• We conjecture that if the Pick matrix II has m negative eigenvalues then there 
exists an interpolating function within a set QV g where g has m poles or m zeroes 
in C_|_. 

• Parameterize all QVE interpolating functions. 

• Characterize the Nevanlinna-Pick interpolation problems solvable by Odd and by 
QVE functions. 

One can then look for generalizations. For example, 

• Formalize the extension of the study of QV g functions, to the cases of: (i) not 
necessarily rational, (ii) matrix valued. 
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• Formalize the extension of the study of the even-odd partitioning of QV functions, 
to the cases of: (i) not necessarily rational, (ii) matrix valued. 
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